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We discuss three questions related to the critical behavior of hierarchical spin models: 1) the hyperscaling 
relations in the broken symmetry phase; 2) the combined use of dual expansions to calculate non-universal 
quantities; 3) the fine-tuning issue in approximately supersymmetric models. 



1. Introduction 

The Monte Carlo method has played a major 
role in the understanding of quantum field theory 
problems for which there are no known reliable 
series expansions. However, the errors associ- 
ated with this method usually decrease like i~^/^ 
where t is the CPU time used for the calculation. 
This often prohibits obtaining high-accuracy re- 
sults. If in the next decades a better knowledge of 
the fundamental laws of physics has to rely more 
and more on precision tests, one should comple- 
ment the Monte Carlo methods with new com- 
putational tools which emphasize numerical ac- 
curacy. This line of reasoning is a motivation to 
use "hierarchical approximations" (an example is 
Wilson's approximate recursion formula ||I|) as a 
starting point since the RG transformation can be 
handled easily. Before attacking the hard prob- 
lem of the improvement of such approximations, 
we would like to show that they allow spectacu- 
lar numerical accuracy, namely errors decreasing 
like e~^*" , for some positive constant A of order 
1 when t is measured in minutes of CPU time and 
0.5 <u< 1. 

For definiteness, we consider Dyson's hierarchi- 
cal model where the total spin in boxes of size 2' 
are coupled with a strength (|)', with c a free pa- 
rameter set to 2^^^"^^!^ in order to approximate 
a nearest neighbor scalar model in D-dimensions 



(see Rcf. |2| for details). These interactions are 
obviously not ultralocal. Models with Fermi fields 
1^ can be constructed similarly by replacing D — 2 
by I? — 1 in c. In addition one needs to specify a 
local measure Wo((?!)), for instance of the Landau- 
Ginzburg type (H/"o(0) = e-^'^'''^'^^) or of the 
Ising type. Under a block spin transformation, 
the local measure changes according to 



W„+i(0)(xe^(4) -* 



(1) 



This recursion formula can be reexpressed in 
Fourier representation as 



i?„+i(fc) « exp(--/3— )(i?„(^))2 



(2) 
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where Rn{k) is the Fourier transform of W^n(0) 
with an appropriate rescaling. It was found 
P] that polynomial approximations of reasonably 
small degree provide very accurate results in the 
symmetric phase. The initial coefficients are the 
only integrals which needs to be calculated, after 
we only perform algebraic manipulations. The ef- 
fects of finite dimensional truncations decay faster 
than e~"^i^* and the finite-size effects like e""*^* 
Direct fits of the susceptibility give a value Q of 
the critical exponent 7 == 1.299140730159. AU 
the digits of this result are confirmed by a cal- 
culation of the eigenvalues of the linearized RG 
transformation about the non-trivial fixed point 
§ofEq. (I). 



Unfortunately, is not possible to get rid of the 
finite size effects in the broken symmetry phase 
iQ] by iterating Eq. (g) enough times because as 
one moves away from the fixed point, rapid os- 
cillations appear. Namely, Rn{k) ~ cos(mc"/^fc), 
where m is approximately the magnetization at 
zero external field. It is easy to see from Eq. 
(0), that for n large enough, the polynomial ap- 
proximation breaks down as the argument of the 
exponential becomes too large. It is nevertheless 
possible to take advantage of the short number of 
iterations for which the low temperature scaling 
is observed to obtain reliable extrapolations, first 
to infinite volume and non-zero external field and 
then to zero external field. Proceeding this way, 
one can calculate the connected q-point functions 
at zero momentum G'q(O), for (7 = 1, 2 and 3 and 
for various values of (3. The scaling and hyper- 
scaling relations together with 77 = imply 



G^(0)oc(/3-/3c)(^/2)(^(i-«/2)-g) 



(3) 



where 7 is the value calculated in the symmet- 
ric phase. Our numerical estimate of the expo- 
nents |q| agree with the predicted values of the 
exponents with an accuracy of 10"'^, which is re- 
spectable when compared to results obtained with 
conventional methods, but large compared to the 
10""'^^ errors obtained in the symmetric phase. In 
order to improve the accuracy of these results, one 
could increase the degree of the polynomial used 
for the approximation, however a short calcula- 
tion shows that with this procedure, the errors 
only decrease like 1/t. A more appealing possi- 
bility consists in factorizing Rnik) into two part: 
one rapidly oscillating and which requires an ex- 
act treatment and another slowly varying which 
can be approximated by polynomials. This sec- 
ond possibility is presently under investigation. 

We now discuss the possibility of replacing the 
numerical evaluation (which has to be repeated 
for each choice of the input parameters) by a suit- 
able series expansion (which could be done once 
for all). Our goal is to obtain an analytical for- 
mula for the non- universal quantities {Aq, Ai . . .) 
appearing in the parametrization of the suscepti- 
bility: 

X^Wc- PyiAo + Ai(/3, -p)^ + ...). (4) 



The linearization method used to calculate the 
universal exponents 7 and A does not provide 
a way to calculate the non-universal quantities 
{Aq, Ai . . .). Indeed, their values "build up" dur- 
ing the crossover between the unstable IR fixed 
point and the high-temperature (HT) fixed point. 
A proper description of this crossover requires 
non-linear expansions about both fixed points. 
This problem can be solved completely for the 
simplified recursion relation for the magnetic sus- 
ceptibility 



Xn+l=Xn + (/9/4)(c/2) 



An ? 



(5) 



which was used in Ref. [pi to estimate the fi- 
nite volume effects. In this simplified model, one 
can construct non-linear functions y and y which 
transform covariantly (multiplicatively) under a 
RG transformation. These two functions are ex- 
pressed as an expansion about the IR and HT 
fixed point respectively. It turns out that there 
exists a duality relation between the series ex- 
pansions of the two functions. It would be very 
interesting to see if similar methods could be ap- 
plied for non-linear sigma models or gauge theo- 
ries. The quantity Aq is a function of (/3c — P) 
with a discrete scale invariance and it can be ex- 
panded in Fourier modes as in the original model 
g. For practical purposes, the contribution of 
the non-zero modes is exponentially suppressed 
and in very good approximation: 



1 



AZa 



^" = h^/ (rf^A)^^y[i-y"'W] 



(6) 



where A is the eigenvalue associated with the un- 
stable direction of the IR fixed point. The lower 
value Za of the integration interval is arbitrary 
and we can choose it at our convenience. We 
have compared the approximate values ^0(771, rh) 
obtained from expansions with m terms for y~^ 
and m terms for y with an accurate value of Aq 
and found that the errors go approximately like 
exp(— _?i'i(r7T, -|- rh) + K2{m — rh)'^). This implies 
that for m + rh fixed, it is very advantageous to 
pick the "self-dual" option m "^ rh. 

We now discuss the fine-tuning question. The 
fact that the bare mass of a scalar field the- 
ory requires a fine-tuning in order to keep the 



renormalized mass small in cut-off units is usu- 
ally regarded as an argument against fundamen- 
tal scalars. A possible resolution of this feature 
consists in adding degrees of freedom in such a 
way that the quantum fluctuations cancel, mak- 
ing small scalar masses a more natural outcome. 
In a recent preprint [^ , we presented two models 
with an approximate supersymmetry and which 
can be solved non-perturbatively with the numer- 
ical methods discussed above. The bosonic part 
of these models is Dyson's hierarchical model. 
The free action for N massless scalar fields (f>x 
reads 



rt/ree _ 






: -^xy'i 



i) 



(7) 



where x and y run over the sites and i from 1 to 
N . The action for free massless fermions reads 



4'^'^^^E^x^^..^^' 



i) 



(8) 



x,y,i 



where the ipx 



Ti^) 



and %px are Grassmann numbers. 
^„ and the 



The explicit form of D^y, its square D^. 
supersymmetry of the free action are given in [p| 
The Grassmann nature of the fermionic fields 
restricts severely the type of interactions allowed. 
For instance, for one flavor the most general local 
measure (without any supersymmetry considera- 
tions) is 



Wi(l),ip,^) ^ w {<!)) + ij4>Ai<i)) 



(9) 



A more interesting local measure with two-flavors 
and a R-symmetry is given in p| . The RG trans- 
formation can be expressed as a finite number of 
convolutions which can be calculated accurately 
in Fourier transform with polynomial approxima- 
tions as in the bosonic case (see ^ for explicit 
formulas). Some numerical results are summa- 
rized in Fig. 1 where the renormalized mass is 
plotted as function of the bare mass (both in cut- 
off units). One sees that in the purely bosonic 
interacting model (LG), the quantum corrections 
to m^ are positive while for a Gaussian model 
coupled to fermions with a Yukawa coupling, the 
corrections are negative. For visual reference, the 
bosonic Gaussian result, where the two quanti- 
ties are obviously equal, is also plotted (straight 
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Figure 1. The renormalized mass {mn) as a func- 
tion of the bare mass (ms) for various models. 



line) . We would like to know if by combining the 
two interactions in an appropriate way, it possi- 
ble to cancel the two corrections. At the one-loop 
level, a simple relation {gy = 8A4) guarantees that 
rrifj goes to zero when m^ goes to zero. How- 
ever, we found numerically that for such a choice, 
m^ ~ 0.044 (in cutoff units) when m^ goes to 
zero. It is possible to fine-tune gy in order to get 
nifi — and the exact critical value of gy is about 
50 percent larger than the perturbative one. We 
are presently trying to build models with better 
naturalness properties. 
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